arXiv:l509.00140vl [astro-ph.HE] 1 Sep 2015 


Extend the Kompaneets Equation to Down-Comptonization Process in a Novel Way 

Xu Zhang" *, Xurong Chen" 

"Institute of Modem Physics, Chinese Academy of Sciences, Lanzhou 730000, P.R. China 
^University of Chinese Academy of Sciences, Beijing 100049, PR. China 
(Dated: September 2, 2015) 

Comptonization is a very important phenomenon in astrophysics. Kompaneets equation describes up- 
Comptonization process in nonrelativistic astrophysics, while it fails for down-Comptonizadon Scattering which 
is the most important radiative transfer process in hard X-rays and y-rays astronomy. In this study, we explore 
both up-Comptonization and down-Comptonization processes. A new relativistic corrections is introduced to the 
Kompaneets equation which is valid in both nonrelativistic energy regime with the photon energy hv « m e c 2 
and the electron temperature KT e « m e c 2 . Numerical comparisons are presented which show excellent agree¬ 
ment between the Kompaneets equation and the new equation. 
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I. INTRODUCTION 


The radiative transfer for photons that results from electron scattering-plays a crucial role in determining the spectra that 
emerge from cosmic X-ray and y-ray sources, and that has been an important topic in both astrophysics and radiation physics. 
With the rapid development of X-ray astronomy in recent years, it seems necessary to improve the available theory for the study 
of the radiative transfer process. 

Comptonization is of vital importance in astro-particle physics. When the high energy photons, usually X-rays or Gamma- 
rays, from the star or some other celestial bodies are going through the media surrounding them, these media would be ionized 
into plasma. Comptonization is just the interaction process between the high energy photons and the plasma. 

In Comptonization, there would exist many Thomson scatterings, many Compton scatterings and many electron-position pair 
creations as well. Since the Thomson scattering cant change the photons energy, the Compton scattering is the majoring process 
in non-relativistic Comptonization. A mono-Compton scattering is just a collision between photo and electron, so we have 
kinetic energy and momentum energy conserved. 

The various studies of the radiative transfer for photons have employed Monte Carlo calculations 11-31 and solution of the 
Kompaneets equation 141 . Monte Carlo calculations is cumbersome and poorly suited for studying the nonlinear problem^ 5-61 , 
e.g. the Bose-Einstein spectrum. With the Kompaneets equation one can handle nonlinear problem. 

The Kompaneets equation is particular form of a Fokker-Plank equation. This equation as follows, 


dn KT e 13 4 dn 

(-Ei)c = - ^N e o- T c-r—{x [— +n(n + 1)]} 

dt m e c 2 x z ox ox 


(1) 


where x= 3v/kTe is the dimensionless photon energy; N,, is the number density of the scattering electron gas; cr T is the Thomson 
cross-section; n(x,t) = n(v,t) is the frequency distribution function of the photon gas. The Kompaneets equation describes the 
up-Comptonization Scattering of low energy photons of frequency v on a dilute distribution of nonrelativistic electrons when 
all photons and electrons are distributed isotropically in their momenta. And the Kompaneets equation is applied with the 
photon energy hv « m e c 2 and the electron temperature hv « KT e . While it fails for use in down-Comptonization of high 
energy photons passing through electron plasma which is the most important radiative transfer process in hard X-rays and y-rays 
astronomy. 

Based on Fokker-Plank equation, Ross and McCray 0 obtained the Ross-McCray equation to describe Compton Scattering: 
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For correct diffusion equation, when photon-gas reaches a thermal equilibrium with the electrons, = 0 should be satisfied. 
While inserting Plank distribution function n(x) = (e x - 1) 1 to Eq.(2), ^ + 0. 

Liu et al. extended the Kompaneets equation [8| which aims at describing a more general Compton scattering process, 
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Following Kompaneets, they assumed that Av is also the small quantity when hv » KT e , and used Av expanding the 
distribution function. While when hv » KT e , assuming Av as the small quantity is not strictly satisfied. The change of the 
photon energy in each collision is given by the following well known formula (if the electron is approximately motionless 
compared with photon before collision, hv » KT e ) 


liAv — hv 


T(1 — cosO) 

1 + /!(1 - cosO ) 



(4) 


The change of the photon energy depends on hv and the scattering angle. With the increasing of the scattering angle, liAv also 
increases. When scattering angle reaches to n , hAv reaches to its maximum value. 

For example: 


• hv =1 KeV, hAv max = 0.0039 KeV; 


• hv =10 KeV, hAv max = 0.38 KeV; 


• hv =30 KeV, hAv mtIX = 3 KeV 

So in the down-Comptonization process KT e « hv, the condition IrAv « KT e is not always satisfied. 

The purpose of this paper is to re-derive the Kompaneets equation in a novel way. We introduce relativistic corrections to 
Kompaneets equation to describe radiative transfer process in hard X-rays and y-rays astronomy. While different from Kompa¬ 
neets, we chose the change of the electron’s momentum \Ap\ expanding the distribution function. In the following this paper, we 
can see A p is more approximately to the small quantity than Av both when hv » KT e and when hv « KT e . 

In Section 2, we will obtain a new relativistic corrections to Kompaneets equation in a novel way. Then in Section 3 we 
compare our new equation to the Kompaneets equation and Liu’s equation. 


II. EXTENTIONS TO THE KOMPANEETS EQUATION 


Following Kompaneets, we treat the radiation as a closed system that consists of photon-gas and electron-gas. Although the 
system can’t be described by a characteristic temperature before the thermal equilibrium being established, the electron-gas itself 
is already in thermal equilibrium as the interaction between electrons is the Coulomb long-range force. For a tenuous radiations 
field, Fermi distribution of electron-gas approximate as Boltzman distribution 


f(p ) = foexp[- 


KT e 2m e 


] 


(5) 


Considering photon which is boson that is in the system, which doesn’t appreciably disturb the field that electron-gas has 
established, but which is capable of absorption and radiation of energy of all frequencies. Over a sufficiently long time, the 
absorption and radiation of the photon-gas caused by the scattering of electron and photon will lead to thermal equilibrium. 

In the electron and photon collisions process that ( p,v,n) —> (//, v', n') leads to a decrease of the photon numbern(v, t), as 
describing the number of electrons in the interval p - p + dp is N e f(p)d 3 p, because the photon is boson, the total transition 
number in unit volume is 


n( 1 + n')N e f(p)d 3 pdw 


( 6 ) 


The inverse process that (//, v', n') —> {p, v, if) leads to a increase of the photon numbern(v, f), the total transition number in 
unit volume is 


n'{ 1 + n)N e f(p')d 3 p'dw 


( 7 ) 
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Where n = n(v, t), n' - n(v', t ) are the photon numbers before and after collision. And p, p' are the electron momentum before 
and after collision. The transition probability dw is same as in collisions process ( p,v,it) —> (/?', v',n') and inverse process that 
(p', V, n ') —> (p, v, it), because in the non-relativistic limit, the Compton differential scattering cross section can be approximately 
expressed by the Thomson section 


• o 2 

dw — cdcTj — c —(1 + cos~9)2nsin0dQ 


( 8 ) 


the Klein-Nishina cross-section ^ has the same value both for the scattering angle 6 and n - 8. Therefore the change of the 
distribution function n(v, 1) that result from photon and electron scattering is 


dn 

(757 )c = ~N e 

at 


fd 3 p‘ f[n(l + n')f(p’) -»'(!+ n)f(p*)]dw 


(9) 


The law of energy and momentum conservation in the non-relativistic approximation are written in the form 
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—n + p — — n' + p' 
c c 

p 2 , p' 2 

hv + -— = hv + -— 

2 m„ 2 m„ 


( 10 ) 


Where it and n' are the direction of photon before and after collision, respectively. The expression Av = v - V , A p — p - p' can 
be obtained from Eq.(10). Retaining only the first order of Av, we obtain 
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( 12 ) 

(13) 


Where y is the angle between p and A/7. For convenience, in the following calculation, we use A p replacing |A/7|. 

By expanding ti — n(v',t) and f(p') in A p to second order, and replacing the frequency v by a convenient dimensionless 


frequency x, where x 


hv 

KT e 


, we obtain 
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The first order is same as expanding n' = n(v\ t) and /(//) in terms of Av. While there are difference for the second and higher 
orders. 

The Compton Scattering and the inverse Compton Scattering consist of the whole scattering process of photon and electron. 
In Compton Scattering, which is satisfied hAv < 0. Using Eq.(12) and (13), we obtain 


1 . hvci 7 . ^ -> hv (hAv) 2 

—2|/7||A/7|cosy = - [- ‘-(n - n')] + (1 - nn')( - ^)hAv + - - T 

m e c z m e c~ 2m e c~ 


2 m,. 


( 16 ) 


The last term is a small quantity compared to other terms, therefore can be ignored. 
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As the second term (1 - nn')(Yy 2 i)hAv < 0, we obtain 

1 


2m, 


2\p\\Ap\cosy\ < | - (ft - n')]\ 


If the electron is approximately motionless compared with photon before collision (hv » KT e ), we obtain 

hvcp ^ (hv) 2 

|-[—5-^(1 -tin')] 

m e c m e c 

| - -— 2\p\\Ap\cosy\ « hAv 
2 m e 

If the thermal energy of electrons is markedly larger than the energy of photons in radiation field! KT e » hv): 

I - [——- n')] hAv 


(17) 

(18) 
(19) 


m e c- 


2m. 


■2\p\\Ap\cosy\ < hAv 


( 20 ) 


So we have proved in Compton Scattering, A p is a smaller quantity than Av in both cases liv » KT e and hv « KT e . 

In inverse Compton Scattering which only happened when low energy photon collides with high energy elcctron(K7',, > hv) 
that is satisfied hAv > 0. 

from Eq.(ll), we obtain \cp(n - n')\ > |(1 - hn')hv\. So that \cp(h - n')\ » |(1 - Ftn')hAv\. Therefore Eq.(16) can be be 
approximately to 


hvcp 


-»— 2\p\\Ap\cosy « -[- j(rt-n')] ~ hAv 


( 21 ) 


2 m e ~ “ m e c 

Eq.(21) is what Kompaneets used when he derived his equations in up-Comptonization process which is applied to non- 
relativistic astrophysics problems. Here we have proved A p is a smaller quantity than Av in both cases hv » KT e and hv « 
KT e . 

So we expand n' = n(V,t) and f(p') in terms of A p to higher order: 
dn , 


n =n + -(--^—^—IpcosyAp + (--J—)Z—(Ap) 2 ] + }-^{[--J—Z—2pcosyAp] 2 
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Insert Eq.(14) and Eq.(15) into the Eq.(9), we obtain 

dn dn 

( dt )c = NA fo +n(n+1)] 
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(24) 


Next, Follow Kompaneets’ method, we first calculate the second integral of the Eq.(24). The other is determined from the 
condition that the equation ought to guarantee conservation of the total number of quanta in the scattering.Let: 


I = 


d 3 p f dWf(p))[hAv + 2—|A/7| 2 ] 2 
J J 2 m e 


(25) 
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Inserting Eq.(13) into the Eq.(25), we obtain 

dWf(p)){(hAv) 2 + 

hv hAv hAv T 9 1 hv 9 , ^ hv hAv hAv , 

-2 (n-n')n'(-) -+(- ) 2 ] 2 + -(hAv)[(—) 2 (n- n'f -2(H -n')n'( —) -+(-) 2 ]} 
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Let: 


h = 

f d3f . 

J " dWf(p))(hAv) 2 

h = 

f di} . 

r 1 hv 9 . -29 
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- 2 (ft - 
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Calculating those integral above, we obtain 


, r 2 KT en x2 ,_ , (hv ) 4 7 n 
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m e c* 

, (/*v ) 4 7 
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Hence we can get 


r 2 KT e „ 28 KT e , . 

/ =- —(hv)~cr T c + — - hvfcr T c 


The Eq.(24) obeys a sort of conservation laws: 


dt c dx 


5 ( m e c 2 ) 2 
d(x 2 j) 


(27) 


Where j is the flow of quanta in the frequency space. Using spherical coordinates (y, 0, pi) to replace (xi, x 2 , X3), Eq.(27) can be 
written as 


A _ 

V r\ )c J r\ 

Ot X ox 


(28) 


The Eq.(24) is of second order relative to x, and dependent on the second derivative '-pi linearly, so that the current must contain 


dx 2 


the first derivative p. Meanwhile in the state of thermal equilibrium, the distribution function is Planckian, n(x) = (e x - 1 ) , 
the flow vanishes. 


dn _ 
dx 

Therefore j follows the form: 


Thus f x = -n(n + 1). 


dn 

j(x) = g(x)[— + n(n + 1 )] 

OX 


(29) 


Inserting Eq.(29) into the Eq.(28), we obtain 


&- 


• S WI^ + 2(„ + iA 

ox z ox 


OX X ox 


(30) 


Meanwhile the Eq.(24) can be written as 
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dn s 1 dn 

—) c + n(n + 1 )]// + 

dt KT e dx 


N ? ( 1 92 r d 2 » 
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+ 2 (n + 1 ) 


dn 

dx 


V 


(31) 



6 


n2 

Comparing Eq.(31) with Eq.(30) and noting that the coefficient of should be the same, g(x) is obtained as: 


g(x) = - = -Ax 2 ( 1 + Bx) 

2 A 1 p 


(32) 


Where A = ^^NeO-jc, B = y yy. Inserting Eq.(32) into the Eq.(29), we obtain Eq.(33), 


(In KT e 


1 d 


14 KT e dn 


(-7T)c=- ? AUr r cy — {x (1 + ^ /L 

of m e c~ x A ox 5 m e c z ox 


x)[ ——h n(n + 1)]) 


(33) 


Eq.(33) is a new equation. Here we should point out that the we use A p rather than Ay in the Taylor expansion. This is the 
main reason that our result is different from the previous results. The new equation can be applied in the nonrelativistic energy 
regime with the photon energy hv « m e c 2 and the electron temperature KT e « m e c 2 to describe a more general Compton 
Scattering process. The comparison between hv and KT e is no longer necessary. When the energy of photons hv is low, then the 
term y ppp x of Eq.(33) can be ignored , so it return to the classical Kompaneets equations. While when the energy of photons 
hv is very high, the effect of the term y jpAx of Eq.(33) could have large effect and hence is not negligible. 


III. NUMERICAL CALCULATIONS 

We apply our new equation to explore evolution of typical X-ray spectra in astro-physics. We calculate Eq.33 numerically by 
using the three-point finite difference method. We fix N e = 2 x 10 ,6 c/?) 3 , which are the typical values for plasma in the accretion 
disk around a compact star. 

The expression of the spectral intensity for an emission line with normal Gaussian profile is given as: 

4/«2 , 

7(v') ~ exp[-——^(v - v 0 ) ] (34) 

(Av) z 

Using the relation /,, = yp-n(v). The initial line profile of a stable source is given as: 

n(x, 0) = f(x) ~ x 3 exp[-—-r(x - * 0 ) 2 ] (35) 

(Ax) z 

In the following calculations, we take the diffusion time scales is T = l.Ox 10 2 . Fig[l]shows our new equation compared with 
the Kompaneets equation in up-Comptonization process of low energy photons passing through electron plasma. We can see the 
difference of the resultant line-profiles between our work and Kompaneets is very small, which confirms that the correctness of 
Eq. (33) is valid in up-Comptonization process. 




FIG. 1: (Color Online) Left :hv 0 = 1.0 KeV, the FWHM of the line as hAv = 0.1 KeV, 0.001 KeV < hv < 5 KeV. Right: hv 0 = 2.0 KeV.the 
FWHM of the line as hAv = 0.2 KeV, 0.001 KeV < hv < 5 KeV. 

For down-Comptonization process, the Kompaneets equation does not work any longer. Here we compare our equation with 
Liu’s equation. Fig. [2] shows our work compared with Liu’s work in down-Comptonization process of high energy photons 
passing through electron plasma. From Fig. [2} we conclude that when the centroid energy as /jvq is not very large, the difference 
between our equation and Liu’s equation is very small, while with the increasing of hv o , the difference is more and more obvious. 
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FIG. 2: (Color Online) Top Left: hv 0 = 10.0 KeV,the FWHM of the line as liAv = 1.0 KeV, 1.0 KeV < hv < 19KeV.; Top Right: hv 0 = 30.0 
KeV,the FWHM of the line as hAv = 1.0 KeV, 11.0 KeV < hv < 49.0KeV; Bottom: hv 0 = 50.0 KeV,the FWHM of the line as hAv = 1.0 KeV, 
11.0 KeV <^v< 79.0 KeV. 


As we know if higher order term of the distribution function is ignored, the small quantity is smaller, the result is more precise. 
When hv o is not very large, the second and higher order term of the distribution function expanded by A p is approximately to 
that of expanded by IrAv. While when hv o is very large, the obvious difference between our equation and Liu’s equation comes 
from the second and higher order term of the distribution function which is expanded by A p is much smaller than that expanded 
by Av. 


IV. CONCLUSIONS 

In this study, we have explored both up-Comptonization and down-Comptonization processes using relativistic corrections 
to Kompaneets equation to describe a more general Compton Scattering process. While different from Kompaneets derived his 
equation, we have choose A p as the small quantity. When hv » KT e , A p is more approximately to the small quantity than Av, 
which ensure the equations is more accurete for down-Comptonization processes. In up-Comptonization process, our results 
match the Kompaneets equation. With introducing relativistic corrections to Kompaneets equation, the non-relativistic process 
with photon energy hv « m e c 1 2 3 4 5 and electron temperature KT e « m e c 2 can be described, with the conditions hv « KT e is not 
necessary. 
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